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f*"*- , We propose a model in which ultra high energy cosmic rays are produced by 

CO 

collisions between neutron stars and axion stars. The acceleration of such a 

m 

cosmic ray is made by the electric field, ~ 10 15 (B/W 12 G) eV cm , which 

o 

is induced in an axion star by the strong magnetic field B > 10 12 G of a 

neutron star. As we have shown previously, similar collisions generate gamma 
P-t. 

<D , 

c~| , ray bursts when the magnetic field is much smaller, e.g. ~ 10 G. If we 
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assume that the axion mass is ~ 10 9 eV, we can explain huge energies of 



the gamma ray bursts ~ 10 54 erg as well as the ultra high energies of the 
cosmic rays ~ 10 20 eV. In addition, it turns out that these axion stars are 
plausible candidates for MACHOs. We point out a possibility of observing 
monochromatic radiations emitted from the axion stars. 
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The origin of ultra high energy cosmic rays ( UHECRs ) is one of most mysterious puzzles 
in astrophysics Q. UHECRs with energies ~ 10 20 eV can not travel in distance more than 
50 Mpc owing to the interactions between UHECRs and the cosmic backgrand radiations. 
Observations, however, show that there are no visible candidates for the generators of such 
UHECRs in the arrival directions of UHECRs. On the other hand, dark matter in the 
Universe |§ is one of most mysterious puzzles in cosmology. Axion |3|,|4j] is a plausible 
candidate for the dark matter. Probably, some of axions may form boson stars ( axion stars 
) in the present Universe by gravitational cooling || or gravitational collapse of axion clumps 
formed at the period of QCD phase transition |3j . Here we explain a generation mechanism 
of UHECRs and discuss a production rate of UHECRs assuming collisions between axions 
stars and neutron stars. We also show that the collisions cause emission of observable 
monochromatic radiations from the axion stars. 

We have previously proposed |7j a possible generation mechanism of gamma ray burst ( 
GRB ). According to the mechanism the collision between axion star and neutron star gen- 
erates a gamma ray burst; the axion star dissipates its mass energy j7||| very rapidly under 
the strong magnetic field of the neutron star. Thus, the energy released in the collision is 
given by the mass M a of the axion star. Typically it is given by M a ~ 10 _5 M Q (10 _5 eV/m a ) 
where m a denotes the axion mass. In the analysis we have taken the mass such as m a ~ 10 -5 
eV as suggested observationally in standard invisible axion models. The mass, M a ~ 10 49 
erg, corresponding to this choice, however, is not enough to explain a huge energy observed 
in GRB980123 even if jet is assumed in the GRB. 

In this paper, we explain an acceleration mechanism of UHECRs assuming the axion 
mass, m a ~ 10~ 9 eV, although the choice is not conventional [@||. The essence of the 
acceleration is that a strong electric field ~ 10 15 (5/10 12 G) eV cm -1 is induced in the axion 
star when it is exposed to the magnetic field B > 10 12 G of a neutron star. This electric 
field can accelerate charged particles and makes them obtain the huge energies ~ 10 20 eV. 
Furthermore, it turns out that the energy ~ 10 54 erg observed in the GRB can be also 
explained assuming a moderate jet of the GRB. This is because the mass of the axion star 
is given by ~ 10 _1 A/ Q ~ 10 53 erg with the assumption of m a . Thus the collisions between 
the axion stars and the neutron stars are possible sources of both UHECRs and GRBs. 



Additionally, it turns out that the axion star is a plausible candidate for MACHO [1C 



because the value of the mass is suitable for explaining the observations of MACHOs. Since 



all of baryonic candidates for MACHOs have been argued to have serious difficulties |TT 



nonbaryonic ones like the axion stars are favored. We discuss that the collisions generate 
monochromatic radiations with a frequency m a /2-K ~ 2.4 x 10 5 Hz. With the detection of 
such radiations we can test our model and determine the axion mass. 

As we shall show below, we need to separate two cases of the collisions, the collision 
with a neutron star possessing relatively weak magnetic field ~ 10 10 G and the one with 
a neutron star possessing relatively strong magnetic field > 10 12 G. In the first case the 
axion star collides directly with the neutron star and dissipates its whole energy inside of 
the outercrust of the neutron star J7|. The ultra high energy cosmic rays are not produced 
in this case but gamma ray bursts are produced. On the other hand, in the second case 
the axion star never collide directly with the neutron star because it evapolates before the 
collision. This is because the axion star induces a stronger electric field |7|;[12| when it is 
exposed to the stronger magnetic field of the neutron star. But such a strong electric field 
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decays very rapidly owing to electron-positron pair creations [pl| when the strength of the 
electric field goes beyond a critical value. Namely, as the axion star approaches the neutron 
star, the strength of the magnetic field around it increases gradually and reaches a critical 
value beyond which the electric field induced in the axion star decays very rapidly. It means 
that the axion star itself decays very rapidly. In this case the ultra high energy cosmic rays 
are produced. Gamma ray bursts might be also produced, but they are emitted in a cone 
with much small solid angle and their duration is very short ( less than millisecond ). 

Let us first explain briefly axion stars ( ASs ) and how they induce strong electric fields 
under a magnetic field of a neutron star. The AS is a coherent object of the real scalar field 



a(x) describing the axion. It is represented by a solution [H,I5| of the equation of the axion 
field coupled with gravity. An approximate form of the solution [|1| is given by 



a {%) = fpQa sm(m a t) exp(-r/R a ) , (1) 

where t ( r ) is time ( radial ) coordinate and Jpq is the decay constant of the axion. The value 
of fpQ is constrained Q conventinally from cosmological and astrophysical considerations 
0J4J] such that 10 10 GeV < fpQ < 10 12 GeV ( the axion mass m a is given in terms of 
fpQ such that m a ~ 10 7 GeV 2 / fpQ ). However, when we assume unconventionally entropy 
productions below the temperature 1 GeV in the early Universe, we may be released from 
the constraints ||16|| . Hereafter we assume that fpQ ~ 10 16 GeV or m a ~ 10~ 9 eV. 

In the formula, R a represents the radius of an AS which has been obtained [§,15] numeri- 
cally in terms of mass M a of the AS; R a = 6Am 2 pl /m 2 a M a ~ 1.6 x 10 5 cmm^" 2 (1O _1 M /M a ), 
with nig = m a /10~ 9 eV and m p i is Planck mass. Similarly, the amplitude a in eq([j]) is 
represented such as a = 1.73 x 10 2 (10 5 cm/R a ) 2 nig . Therefore, we find that the solution 
is parameterized by one free parameter, either one of the mass M a or the radius R a . It is 
also important to note that the solution is not static but oscillating with the frequency of 
m a /2ir. It has been demonstrated JT7| that there is no static regular solution of the real 



scalar massless field coupled with gravity. This may be general property of the real scalar 
massive field. 

The AS mass is determined by physical conditions under which the AS has been formed; 
how large amount of cloud of axions are cooled gravitationally to form the AS, etc.. The 
situation is similar to other stars such as neutron stars or white dwarfs. A typical mass 



scale in these cases is the critical mass ||18|| ; stars with masses larger than the critical mass 
collapse gravitationally into more compact ones or black holes. In the case of the AS, there 
also exists a critical mass M c which is given by |14] M c ~ 1O~ 1 M nig 1 where M represents 



the solar mass. Therefore, we adopt this critical mass M c as a typical mass scale of the ASs. 
( The critical mass is the maximal mass which ASs can take. Thus, their masses, in general, 
are smaller than this one. Since energies released in GRBs by the ASs are given by masses 
themselves, the maximal energy in GRBs is given by the critical mass ). 

Although the gravitational cooling for star formation is in general ineffective because it 
is too slow process, it has been shown H that the cooling is very effective for the real scalar 
axion field. Thus, the axion stars can be easily formed gravitationally in a gas of the axions. 
It is reasonable to assume that the most of the axions in the Universe forms the axion stars. 

Let us explain how an AS induces an electric field under a magnetic field B of a neutron 
star. Owing to the interaction between the axion and the electromagnetic field described by 
L a77 = caa(x)E ■ Bj fpQir, where the value of c ( of the order of unity ) depends on axion 
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models |T^,^,|], the Gauss law is modified [2TJ such that dE = —cad ■ (a(x)B)/ fpQ7t + 
"matter". The last term "matter" denotes electric charges of ordinary matters. The first 
term in the right hand side represents an electric charge made of the axion. ( This charge 
is oscillating and so there exist a corresponding oscillating current, J a = cad t a(x)B / fpQTi. 
Thus, radiations are emitted by the AS, which might be observable. We will discuss it 
later. ) Thus, we find that the electric field E a is induced; E a = —caa(x)B/fpQ7i with 
a ~ 1/137, Numerically, its strength is given by 

E a ~ 10 15 eV cm" 1 B 12 m 9 , (2) 

with B12 = B/10 12 G, where we have used the solution in eq([I|) for the critical mass. 
Obviously, the spatial extention of the field is given by the radius R a ~ 1.6 x 10 5 m g ' cm of 
the AS. " 

This electric field is oscillating with the frequency, m a /2-K ~ 2.4 x 10 5 mgHz. Thus a 
particle with electric charge Ze can be accelerated in a direction within the half of the period, 
Tr/m a or in a distance ~ R a ( ~ ir/m a x light velocity ) by this field, unless it collides with 
other particles within the period. Thus, the energy AE obtained by the particle is given by 

AE = ZeE a x 7r/m a x light velocity ~ 10 20 Z eV£ 12 . (3) 

Therefore, the electric field can accelerate the charged particle so that its energy reaches 
~ 10 20 Z eV. These charged particles may be baryons or electron-positron pairs produced 
by the decay of the electric field itself as discussed soon below. 

Comment is in order. It seems apparently from eq([J) that stronger magnetic fields 
yield cosmic rays with higher energies. Stronger magnetic fields, however, induce stronger 
electric fields, which are unstable against electron-positron pair creations. Therefore, strong 
magnetic fields do not necessarily yield cosmic rays with higher energies. 

If the particles collide with other particles on the way of acceleration, in other words, 
their mean free paths are shorter than R a , they can not obtain such high energies. It is easy 
to see, however, that the mean free paths of quarks or leptons with much higher energies 
than their masses are longer than R a in magnetosphere of neutron star. This is because 
since cross sections, a, of quarks or leptons with such energies E behaves such as a ~ 1/E 2 , 
mean free paths ~ 1/na is longer than R a ~ 10 5 cm unless number density n of particles 
around the AB is extremely large ( i.e. n > 10 40 /cm 3 ). 

As is well known, the strong electric field is unstable |T3| against electron-positron pair 



creations. This implies that AS decays into the pairs when the electric field induced in AS 
is sufficiently strong. 

Let us estimate the decay rate of the field and show that the AS decays before colliding 
with a neutron star whose magnetic field at surface is stronger than 10 12 G. We also show 
that the AS can collide directly with a neutron star whose magnetic field is relatively weak 
~ 10 10 G. 

The decay rate Rd of the field per unit volume and per unit time is given by 



aE 2 a ™ exp (-m 2 e Tm/eE a ) 

R d=—2-l^ "a ( 4 ) 



7T " rr 



where m e denotes electron mass. The rate is very small for an electric field much weaker 
than m 2 7r ~4x 10 16 eV/cm. The electric field of AS, however, can be very strong and it 
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can be comparable to m 2 7r. Therefore, the rate is much large. Numerically, it is given by 
R d ~ 10 47 B\ 2 m\ cm~ 3 s _1 S^° =1 exp(— 0.7 x 10 2 n/B 12 m g )/n 2 . Since the spatial extention 
of the field is approximately given by 10 5 rrig 1 cm, the total decay rate W of the field in AS 
is ~ 10 62 B\ 2 nig 1 s~ 1 E^iexp(-0.7 x 10 2 n/ B 12 m 9 )/n 2 . Numerically, it reads 

W ~ 10/s for B l2 = 0.5, W ~ 10 6 /s for B l2 = 0.55, W ~ 10 32 /s for B l2 = 1 (5) 

with m 9 = 1. 

Therefore, we find that the AS decays very rapidly ( or almost suddenly ) when it 
approaches a region where the strength of the magnetic field reaches a critical value of about 
10 12 G. Hence, the AS evapolates before colliding with the neutron star whose magnetic field 
at the surface is stronger than 10 12 G. The whole energy of the AS is transmitted to the 
charged particles, each of which can obtain energies ~ 10 20 eV. These particles are emitted 
into a cone with very small solid angle. They form an extremely short pulse whose width 
being less than millisecond. Actually, when we suppose that the relative velocity of the AS 
is equal to light velocity, it decays approximately within a period of 10~ 4 sec ~ 10~ 5 sec; in 
the period it passes the region where the magnetic field increases from 0.5 x 10 12 G to 10 12 
G. These leptons may be converted into baryons and photons through the interactions with 
themselves, interstellar medium or ejection of progenitor of the neutron star. Thus when 
the magnetic field is sufficiently strong, ultra high energy cosmic rays can be produced. 

We comment that the velocity of AS trapped gravitationally to a neutron star, is ap- 
proximately given by the light velocity just when it collides with the neutron star. This 
is because an AS can be trapped to a neutron star when the AS approaches it within a 
distance ~ 10 11 cm as we will show below, and then, the AS collides with the neutron star 
losing its potential energy and angular momentum by emitting gravitational waves. Thus, 
the velocity of the AS reaches approximately the light velocity when it collides with the 
neutron star. 

Furthermore, we can see from eqs(|2|) and (^) that the critical electric field depends on 
the factor of Bi 2 x m 9 . On the other hand, the energy AE obtained by accelerated charged 
particles depends only on the factor of B\2- Therefore, we find that if m a is smaller than 
10 -9 eV, the maximal energy of cosmic rays can be larger than 10 20 eV when B > 10 12 G. 

It is easy to see that the decay rate is negligibly small for the case of weak magnetic field 
~ 10 10 G. Thus the AS collides directly with such a neutron star and dissipates the whole 
energy in the outercrust of the neutron star. Actually, the rate of the energy dissipation 
in the crust has been estimated previously |7|,|8| and given approximately by 10 46 erg/s cm 3 , 
while the energy density of the AS is given by 10 38 erg/cm 3 . Thus even if the velocity of the 
AS is equal to the light velocity, it dissipates its whole energy in the crust: it never reaches 
the core of the neutron star. This estimation has been performed by noting that the energy 
dissipation of the AS arises due to the dissipation of an electric current ( = a c E a ) induced 
in the crust with conductivity a c ; the value of 10 26 /s has been used for <j a p2 . 



Anyway, this collision generates gamma ray bursts with energies ~ 10 53 erg. The ejection 
could be emitted into a cone with small solid angle as a jet. This is because particles ( mainly 
irons ) of the neutron star are accelerated and emitted in the direction of the strong electric 
field ~ 10 13 -Bio eV/cm. The fact that the whole energy of AB is dissipated only in the 
outercrust, implies that the ejection may be particles composing the crust. Thus a fraction 
of the baryon contamination in the jet is less than 10 _5 M Q as required observationally. 
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In the above case, the AS dissipates its whole energy in the first collision. On the other 
hand, an AS may collide several times with a neutron star when its mass is much smaller 
than the critical mass ~ 1O -1 M . The mass has been assumed as a typical mass of ABs. 
We see from the general formula of R a that the radius of the AS with smaller mass than 
the typical one is larger than the typical radius ~ 10 5 cm. For example, if its mass is given 
by 1O~ 2 M , the radius is about 10 6 cm. This is comparable to the radius of the neutron 
star. In addition, an electric field induced in such an AS is weaker than the typical one 
eq(0). This implies that the rate of the energy dissipation in a neutron star is smaller than 
the one quoted above. In such a case the collisions might occur several times. We expect 
that these collisions leads to GRBs with pulses of longer duration and softer gamma rays. 
Observationally the only energies of GRBs with long duration have been measured. We 
predict that the energies of GRBs with short duration are much larger than those of GRBs 
with longer duration. 

We now wish to estimate an energy release rate per unit volume and per unit time; the 
energy released as ultra high energy cosmic rays. In order to do so we assume that the 
dark matter is composed mainly of axion stars. The estimation, however, involves several 
ambiguities associated with number density of neutron stars, energy density of dark matter 
or velocity of ASs in the Universe etc. Therefore the estimation does not lead to a conclusive 
result although our result is consistent with the observations [[!]. 

First we note that luminosity density around our galaxy is observationally given by 
2 x lO 8 hL /Mpc 3 where h is Hubble constant with the unit of 100 km/s Mpc and L Q 
represents solar luminosity. We take a mass density p corresponding to this luminosity 
density such as p ~ lO 9 M h 2 /Mpc 3 . Then, we suppose that the number density of neutron 
stars is given by ~ 10~ 2 p. This comes from the fact that the rate of appearance of supernovae 
is about 0.1 ~ 1 per 10 years in our galaxy and the age of the galaxy is about 10 10 years. 
Probably, the rate of the appearance could be larger in early stage of our galaxy than the 
one at present. Thus, in our galaxy 10 9 neutron stars might be involved corresponding to 
the number 10 11 of stars in the galaxy. Hence we guess that the number of the neutron 
stars are equal to about 10~ 2 times the number of the stars. All of these neutron stars are 
assumed to possess strong magnetic field > 10 12 G. Furthermore, to estimate the rate of 
the energy release we need to know the average density e of the dark matter. Here we use 
a value [[| of 0.5 x 10 _24 g/cm 3 , which represents a local density of our halo. Using these 
parameters we can estimate the rate owing to the collision between the AS and the neutron 
star if the cross section of the collision is found. The collision cross section for a neutron star 
to trap an AS is estimated in the following. Namely an AS is trapped by the neutron star 
when the AS approaches the neutrons star within a distance L c in which its kinetic energy 
M a v 2 /2 is equal to its potential energy 1.5 x M & M a G/ L c around the neutron star. G is 
gravitational constant. Here the mass of the neutron star and the relative velocity v are 
assumed to be 1.5 x M and 3 x 10 7 cm/s, respectively. Thus, the cross section is found such 
as L 2 7r with L c ~ 6 x 10 11 cm. After being trapped, the AS collides with the neutron star 
by losing its potential energy and angular momentum owing to the emission of gravitational 
and electromagnetic waves; electromagnetic radiations arise due to the oscillating current 



J a . The merger is similar to neutron star - neutron star merger |23|]. Time scale from birth 
to merger is much less than the age of the Universe when the distance L c between two stars 
at the birth is given such that L c ~ 10 11 cm. 
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Therefore, we obtain the rate of the collision per Mpc 3 and per year, 



e x l(r 2 p x L 2 c n xv xl year/lO _1 M ~ 3 x 10~ 9 h 2 /Mpc 3 year. (6) 

Since the energy of ~ 10 53 erg is released in the collision, we find that the rate of the energy 
release is given by 

~ 3 x 10 44 h 2 erg/Mpc 3 year. (7) 

which agrees well with the observed one . Taking account of the fact, however, that there 
are several ambiguities in the parameters used above and in the observations, we understand 
that our model can explains roughly the observations M. 

Finally we discuss two possibilities of the observation of the axion stars. One is associated 
with gravitational lensing and the other one associated with monochromatic radiations from 
the ASs. 

Since we assume that the halo of our galaxy is composed of ASs and their mass is ~ 
1O _1 M , the ASs are plausible candidates for gravitational microlenses, i.e. MACHO. Since 
baryonic candidates like white dwarfs, neutron stars, etc. have serious problems, nonbaryonic 



candidates are favored |TT[ . The problems are associated with chemical abundance of carbon 
and nitrogen in the Universe: If these baryonic stars are MACHOs, an overabundance of 
the elements is produced far in excess of what is observed in our galaxy. Hence, the ASs are 
theoretically the most fascinating candidates for MACHOs since they are also candidates 
for the generators of both UHECRs and GRBs. If the most favorable mass of the MACHO 
is O.5M , we need to choose m a ~ 0.2 x 10~ 9 eV since the mass of the ASs is given by 
~ 10~ 1 M & /mg. Smaller axion mass leads to stronger electric field as well as larger mass of 
AS. Thus, it yields higher energies of the cosmic rays, ~ 5x 10 20 eV and of GBRs, ~ 5x 10 53 
erg than the ones we have claimed above. Accordingly, the determination of the mass of 
MACHOs gives the upper limit of both the energies of the ultra high energy cosmic rays 
and the energies released in the gamma ray bursts in our mechanism. 

We point out another way of the observation of the axion stars. Since the electric field 
induced in ASs is oscillating, electromagnetic radiations are emitted ]15 , |24] from correspond- 



ing oscillating electric current J a as mentioned above; the frequence of these monochromatic 
radiations is ~ 2.4 x 10 5 m 9 Hz. We expect that such radiations can be detected in advance of 
UHECRs; they are emitted by the AS revolving around a neutron star before the AS decay- 
ing into charged particles. It is easy to estimate their luminocity |24| ~ 6.7 x 10 41 B\ 2 erg/s. 



If the emission arises in a distance ~ 10 Mpc from the earth, we obtain the flux at the earth, 
~ 10 9 Jy Bf 2 /mg] we have assumed that the velocity of the AS revolving is ~ O.lx light 
velocity. Although the possibility of observing the radiations is very intriguing, it might be 
difficult to detect the radiations with such a low frequency because they could be absorbed 
by interstellar ionized gases. 
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